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Abstract. Recently we proposed a novel approach to the formulation of relativistic dissipative hydrody- 
namics by extending the so-called matching conditions in the Eckart frame [Phys. Rev. C 85, (2012) 14906]. 
We extend this formalism further to the arbitrary Lorentz frame. We discuss the stability and causality of 
solutions of fluid equations which are obtained by applying this formulation to the Landau frame, which 
is more relevant to treat the fluid produced in ultra-relativistic heavy-ion collisions. We derive equations 
of motion for a relativistic dissipative fluid with zero baryon chemical potential and show that linearized 
equations obtained from them are stable against small perturbations. It is found that conditions for a fluid 
to be stable against inflnitesimal perturbations are equivalent to imposing restrictions that the sound wave, 
Cs, propagating in the fluid, must not exceed the speed of light c, i.e., Ca < c. This conclusion is equivalent 
to that obtained in the previous paper using the Eckart frame [Phys. Rev. C 85, (2012) 14906]. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 24.10.Nz Hydrodynamic models 



1 Introduction 

It is known that one faces problems of instability and 
violation of causality in solutions obtained from the rel- 
ativistic Naiver-Stokes (NS) equations (observed in the 
first-order theories [UH] which naively extend the non- 
relativistic NS equation). Israel and Stewart (IS) provide 
a phenomenological framework for a relativistic dissipative 
fluid [2111] accounting for these problems. In their model, a 
possible general form of the non-equilibrium entropy cur- 
rent is described by the dissipative part of the energy- 
momentum tensor and by the particle current up to second 
order in the deviation from the equilibrium state. After the 
stability and causality of the IS theory have been shown 
by Hiscock and Lindblom [SIISII^, the causal dissipative 
hydrodynamical model was adapted to study the dynam- 
ics of hot matter produced in ultra-relativistic heavy-ion 
collisions by Muronga [Hlin]- In recent years, the IS model 
plays an important role in the analysis of the experimental 
data obtained by RHIC and LHC (see, for example, [TU]). 

In parallel to the application of IS theory, investiga- 
tions of the basis of relativistic dissipative hydrodynamic 
theory was continued tlLJl^iS, 14,,15,16,Jii8,I9,. This 
is because the IS theory, in its present form, is too general 
and complex from the point of view of quantum chromo- 
dynamics (QCD), which is believed to be the fundamental 
theory for strongly interacting systems [11, . Another rea- 
son is that the theory of relativistic dissipative hydrody- 
namics is not yet fully understood because, for example, 
the equation of motion of the fluid it uses depends on the 



choice of the Lorentz frame [16] (or on the definition of the 
hydrodynamical flow). Since the dissipative part of the 
energy-momentum tensor ST'^'^ and the particle current 
SN^ cannot be determined uniquely by the second law of 
thermodynamics, one usually introduces some constraints 
to fix them. These constraints are known as matching con- 
ditions. 

The other reason to introduce these matching (or fit- 
ting) conditions [3] is the necessity of matching the en- 
ergy density and baryonic charge density, (e, n), in a non- 
equilibrium state to the corresponding equilibrium densi- 



ties: (e 



eq 5 ' ''cq 



-oqi 
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or, equivalently. 



0. 



(1) 



Matching conditions allow to determine the thermody- 
namical pressure, Peq(£cq, 'T'cq), (defined as work done in 
isentropic expansion) via the equation of state for the equi- 
librium state. Here P^q should be distinguished from the 
bulk viscous contribution, 77 = —^Af^^ST'^'^, present in 
the energy- momentum tensor [4]. Finally, matching con- 
ditions are also needed because they are necessary for the 
thermodynamical stability of the entropy current (see Ap- 
pendix A in Ref.[2D]). However, the matching conditions 
given by eq.(IT]) are not unique. So far, except of some re- 
cent works [roill7| . they were not investigated in detail. 

A state of a relativistic dissipative fiuid is described by 
the energy- momentum tensor, T^'^{x), and by the baryon 
number current, N^^ , which obey the conservation laws 



0, 



(2) 
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K = 0, (3) 

and the second law of the thermodynamics 

Si", > 0. (4) 

Because of the uncertainty in definition of the flow ve- 
locity u^{x) for a non-equilibrium fluid, one needs, unlike 
in the case of perfect fluid, to fix the frame for the fiuid 
considered. Two special frames can be defined: Landau 
and Eckart. 

The Landau frame [2] is defined by the vanishing of 
the energy flow, which consists of heat flow q'^ and net 
baryon number flow V^] 



0, 



P.. 



V^, (5) 



''eq 



with Z\(^ = g^^—u^Ui, being the projection operator orthog- 
onal to the four vector u^. In the Eckart frame, the hydro- 
dynamic flow velocity (with normalization u'^u^ = 1) is 
defined by using baryon charge current = N^/ \/N^N^. 
In this frame one always has 



= A^N^ ^ 0, 



(6) 



Equations of motion for the fiuid should not depend on 
the choice of Lorentz frame. Since the relativistic dissi- 
pative fiuid dynamics with extended matching conditions 
has been already formulated in the Eckart frame [3T] , it is 
interesting to check it in another Lorentz frame, for exam- 
ple in the Landau frame. Also, the stability and causality 
conditions for a relativistic dissipative fluid should not 
depend on the Lorentz frame used. The purpose of this 
article is therefore to investigate how the stability and 
causality conditions should be imposed on a fluid depend- 
ing on the Lorentz frame used. 

This paper is organized as follows: In Sec HI we re- 
formulate relativistic dissipative hydrodynamics with ex- 
tended matching condition (as originally introduced in 
ref. |21] ) in an arbitrary Lorentz frame with ^ and 
7^ (neither the Eckart frame with — nor the 
Landau frame with — 0). In such a general frame, 
the flow velocity field u^(x) may be determined by us- 
ing coexisting and V^. We then choose the Landau 
frame and consider a fiuid with zero baryon chemical po- 
tential lib — 0, appearing in the central rapidity region 
of the ultra-relativistic heavy-ion collisions. In Sec|31 we 
check the stability and causality of the fiuid obtained from 
our model in this frame. We close with Sec|4] containing a 
summary and some further discussion. 



where a = Hb/T, |3^' = with (3 = 1/T, T and 
are, respectively, the temperature and baryon chemical 
potential. In the local equilibrium case it is given by 



(8) 

(9) 



where Tj!^ and N^^ are equilibrium energy-momentum 
tensor and baryon charge current: 



(10) 
(11) 



In this case, the energy-momentum conservation, T^q.^ — 
0, and the baryon number conservation, Nj^^.^ = 0, to- 
gether with thermodynamic relations result in the locally 
conserved entropy current: 



^cq;/^ 



.^A^e^q + [3x;f.T,\^ + [/J^Pcq];p = 0. (12) 



To this current we now introduce dissipative corrections 
by adding corresponding dissipative corrections ST'^'^ and 
SN'^ to the energy-momentum tensor and to the baryon 
number current appearing in 0q and ^^q". In this way one 
extends the expression of equilibrium entropy current to- 
wards the off-equilibrium entropy current: 



X " 



(13) 



= SN^', S^^" = ST^" + Y i—Ac.pST"'^. (14) 



Notice that term proportional to AapST"^ , which appears 
due to the natural extension of eq.®, results in 



(15) 



where 77 is bulk pressure, cf. Eq. ([18)) below. It means 
than that entropy density is not maximal in spite of equi- 
librium state used, the entropy current ea. ([T^ is thermo- 
dynamically unstable [20]. Therefore the term A^pST"-^ is 
usually dropped (i.e., x should be put equal zero, x = 0). 
However, the problem of thermodynamic instability can 
be avoided by simultaneously demanding the natural ex- 
tension of the form of the entropy current feq. (|13p with 
(fH)) ) and the following general extension of matching con- 
ditions; 



X^O, u^uJT^" ^ 0, u^Sm ^ 0. 



(16) 



2 Extended matching condition in arbitrary 
frame 

2.1 General form of ofF-equilibrium entropy current 

The general off-equilibrium entropy current can be written 
in the following simple form using the vector cj)^ and 2- 
rank symmetric tensor <?^'', 



(7) 



2.2 Extended matching conditions 

To restore thermodynamical stability as discussed above, 
we propose to impose the following extended matching 
conditions on the dissipative correction of the energy mo- 
mentum tensor and baryon charge current, ST^^'^ and SN^, 
respectively: 



A, SN^'ua = Sn. 



(17) 
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With these conditions, off-equihbrium contributions for 
the energy momentum tensor and baryon charge vector in 
general Lorentz frame are 



where 6 is the divergence of the flow velocity field, 9 = u?^. 
Using the definition that a — one can also write the 
entropy production as 



ST^"" = Au'^u" - nA^"" + W^v!" + Wu^ + 7^^^ (18) 
= 5nu^ + V^Al, (19) 

where 77 is bulk pressure, tt'^" is shear tensor and 
is net flow of the baryonic charge. In this case, the off- 
equilibrium entropy current eq. (jl3|) is given by 



■-cq 



yl + x7T]-u^. 



(20) 



To ensure thermodynamic stability, one may impose a con- 
dition on the entropy current ^ ea. ipOl) . demanding that 



0. 



(21) 



This requirement can be satisfied by the following unique 
condition. 



X77 = -ii}^5n + A. 



(22) 



dUf^ 

^ dT 

ne + \ 



T 

J dSn dA^ 
dr dr 



V T 

^1 



(27) 



Note that the above equation can also be expressed in the 
following form: 



T 



^ dr 

dr fib 

dSn dA^ 



dr dr 



[fibSn - A]e, 



(28) 



Note that, when both A and Sn are set equal to zero, x 
should also be zero, as so far considered widely in the lit- 
erature. However, once one assumes that A and/or 
5n ^ 0, there exists a term proportional to x- By the ex- 
tended thermodynamical stability condition, eg. ip^ . the 
entropy current in the arbitrary Lorentz frame is 



(23) 



One can obtain the entropy current corresponding to Eckart 
or Landau's formulation in the limit and — ^ 0, 

respectively. The entropy current eq. (|23p can be rewritten 
using eqs. dT5| . ([T^ and thermodynamical stability con- 
dition eg. ((22)) . 



(24) 



Thus the entropy production in this case is given by 

It should be noted here that ea. (j25l) is exactly the same 
as eq.(19) in Ref. [H] obtained in the Eckart frame when 
using extended matching condition. This is because S"/^ is 
a scaler and so it does not depend on Lorentz frame. This 
can be seen in the form of the entropy current, eq. (|24p . The 
expression for has the same form both in Eckart frame 
(VM = but 5N^' ^ 0) and in Landau frame {W^' = but 
PxST^^" ^ 0). 

The entropy production is explicitly given by 



/3 77(? + /3 V^^tiA)^^^ 
. d5n „ dA 



where W^' = W^' - fibV" and V = fibV". 



2.3 Constitutive equations for a dissipative fluid in the 
Landau frame 



igWe shall now consider the Landau frame {W^ = 0) which 
is more relevant in the context of ultra-relativistic heavy- 
ion collisions. In particular, we consider a fluid in the cen- 
tral rapidity region where it is expected that the baryon 
chemical potential fib is small. In this paper, we assume 
that, for simplicity, fib — Q (this implies that Ucq = 0.) In 
this case, the entropy production takes the simple form; 



-ne + v^^ux^t:^'' 
-no + v^^ux)^^^^ 



Axn) 



dr 



{xn)e]. (29) 



In the second line of eq. (pUj) we have used ea. ip^ with 
fib — 0. For the x iii ea.(P^ we then use 



77 



(30) 



i.e., the form of x used in eq.(21) in the ref. [21] with = 
0. The second law of thermodynamics is guaranteed (with 
C, and rf being, respectively, the bulk pressure and shear 
viscosity, which are all positive constants) if the entropy 
production is given by 



772 
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This requirement determines the following constitutive equa- 
tion for, respectively, bulk and shear pressure: 



n , „ , K dn 



(32) 
(33) 



Because eg. ([5^ includes term proportional to 1 /U dU/ dr, 
one can introduce into the bulk pressure 77 an arbitrary 
constant z (with dimension [GeV]"*) and write 

7J d n kC d n c, 

z dr z z clT z z 

We shall now consider small perturbations of 77 fields. The 
bulk pressure 77 can be written as 



77 = 77o + (577 



(35) 



with the background reference field 77o and its perturba- 
tion field 511. One can also regard 77o as the value of 77 at 
initial proper time tq, 77o = 77 (tq). Correspondingly, one 
can write 



obtained from the initial flow vector 



field Uq at tq. In this case, the perturbation field 611 can 
be interpreted as 611 = 77(t) — 77(ro). In this sense, 77o is 
a kind of parameter showing the degree of non-equilibrium 
at initial stage. Identifying the arbitrary constant z with 
77o and noticing that ^ln(I -I- ^) w JT^'El^^n, one can 
rewrite the above equation as: 



Tn- 



dsn 



6n = -ai + x)S0, 



where the relaxation time Tn is given by 
rn - C(2C + «/77o) 



(36) 
(37) 

(38) 



and 69 ~ 9 ~ 9q. It is interesting to note that 77o con- 
tributes to the relaxation time Tn . This means that relax- 
ation processes may depend on the initial condition. Note 
also that if k = then contribution k/ Uq in the relaxation 
time would disappear. 

A similar approach can also be applied to the tt^'^ field, 
with perturbation of the shear viscosity around t:'^'^ — 0, 
leading to 



d6^^^''' 
dr 



+ 6tt'"' = 2r;V 



where the relaxation time r^- is given by 



(39) 



(40) 



3 Stability and causality of the fluid 

The stability of a general class of dissipative relativistic 
fluid theories was investigated by Hiscock and Lindblom 
[51[S1[7]. Denoting by 6V{x) the difference between the ac- 
tual non-equilibrium value of a field V{x) and the value 



■in the background reference state, Vb(a;), we assume that 
variations 6V are small enough so that their evolution is 
adequately described by the linearized equations of motion 
describing the background state. We shall now investigate 
the stability of the fluid obtained in our model following 
a prescription proposed in Ref. [SUHIIZ]- In what follows: 

1. The background reference state is assumed to be ho- 
mogeneous in space. Notice that, unlike in Ref. [illSlIT], 
in our case it is not an equilibrium state but rather a 
non-equilibrium one with 77 = 77o and with tt'"^ = 0. 
Furthermore, the background space-time is assumed to 
be flat Minkowski space, so that all background field 
variables have vanishing gradients. 

2. We consider following plane wave form of perturbation 
propagating in x direction 

6V = 6Vq exp{ikx + Ft). (41) 

Linearized equations for dissipative fluid dynamical model 
are given by 

5[T^'''],,, = 0, (42) 

with the perturbed energy-momentum tensor: 

5[T^'''] = (e*q + p:^)[6u^'u'' + u^'5u'') + (5el^ + K6n)u^'u'' 

- (JP^q + 6n)At"' + 5t:'"' . (43) 

Here e*^ and P*^ are energy density and pressure in the 
background non-equilibrium state 



k77o, P* = Peq + 77o. 



However, since (5[77o] = (it has vanishing gradient and is 
constant in r) , terms proportional to 77o do not contribute 
to the linearized equations eg. (1431) (we ignore terms like 
77o(5m'' in the linearized equation). Hence, one can replace 
in the eq.(|33]) e* and P* by the, respectively, £oq and 



Poq- For baryon charge current, since we deal with a fluid 
in the region where the baryon chemical potential can be 
considered /ib = and the net baryon density ricq — 0, 
one has 5[N^'] = 0. 

The perturbed fluid dynamical fields must satisfy con- 
straints 



Hence, in the rest frame of fluid. 



(1,0,0,0), the 



proper time r component of the flow velocity field van- 
ishes, 6u^ = 0. We therefore obtain the following lin- 
earized equations for the energy-momentum tensor and 
baryon number current: 

^[r'"^];^ = (eeq + Peq)((*fc'5u")^^" + PSun 

+ (Pfeoq + Kr6n)u'' - (VJPcq + ^(577) 

+ {ik)6TT^'' = 0. (44) 

The linearized constitutive equations for 77 and tt'^'^ have 
the form (with k = I + k ) : 

(1 + TnP) 



C 

(1 + T^P) 

2r] 



-5n = —k{ik)6u^, 



(45) 



(^vr^^^ = ~\iik)[6^6u'' + 6^6u'' - h^"'6u^] 



(46) 
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The parameters ^ and ^" introduced above have been ab- 
sorbed in the expressions for relaxation time, tu and r^, 
respectively. On the other hand, the parameter k in the 
expression of the entropy production is kept and not ab- 
sorbed in Tn- Its role will be discussed later. 

All perturbation equations can be expressed in concise 
matrix form: 



M^5Y^ = 0, 



(47) 



where 5Y^ represents the list of fields. The system matrix 
can be expressed in a block-diagonal form when one 
chooses the following set of perturbation variables [B] 

5uy,5TT''y, Su'',Sn-'\6Try'',STTyy -Stt''}. (48) 



In this case, 



Q 



M : 



R 



R 



(49) 



where the matrices Q and R are given by 



Q 



R 



ik 



V 



r 

9£cq 






ikhcq 






rhcq 


ik 


ik 


ikk 


l+TnT 

C 






ik 


i+r„r 

4,7/3 / 



ik 



ik 



(50) 



respectively, and I is the 2x2 unit matrix. The /icq denotes 



the enthalpy density which is defined by h^q — £cq 
For r and k satisfying dispersion relation 

[detM] = [det R]2[det Q] = 0, 



cq- 



(51) 



one has plane- wave solution such as ea. (|4T|) for the lin- 
earized equations of the system ea. (|T7)) . 

In what follows we shall discuss in detail the stability 
of transverse and longitudinal modes separately. 



3.1 Propagation of the transverse mode 

The dispersion relation obtained by setting 

77 det(R) = {T^Kq)r^ + h^qF + T]k^ = (52) 

corresponds to the solution of the perturbation equation 
which is referred to as the so-called transverse mode. The 
solution of the above equation ([5^ is given by 



r : 



/icq ± J hi - Ar]{T^heq)k'^ 



2{T„hcq) 



(53) 



Note that we have always Re[/^] < independent of the 
value of k, which means that any small perturbation prop- 
agating in the transverse direction (perpendicular to the 



X axis, direction which the perturbation wave propagates) 
will be damped with time r. Since the general solution 
is given by a linear combination of those solutions of the 
transverse mode, one can say that the plane- wave solution 
of the mode is stable against small perturbation. 
Note also that when wave number k > kc, where 



fee 



(54) 



the linear perturbation wave propagates towards the trans- 
verse direction, but waves with wave number k < kc are 
damped. 



3.2 Propagation of the longitudinal mode 

Frequencies of the so-called longitudinal mode (propagat- 
ing parallel to the x direction) are given by the roots of 
the following dispersion relation: 



n=4 

det(Q) = 9"^" 

n=0 



0, 



where the coefficients g„ are given by 

94 = ho. 



93 



92 



91 



hon I — + — 



1 
Tn 
1 

TnTr, 



h k(1 - K——^) 

Tn 



de. 



cq 



90 



^h 



3t^ 

1 1 
h - 

Tn Tt, 

ii] 1 
3t^ Tn 

dPcq 1 1 



dec 



OR 



;.2 



9e, 



cq 



cq 



Tn Tn 



^- k 



cq 



dScn Tn T^ 



(55) 

(56) 
(57) 

(58) 

(59) 
(60) 



When all coefficients 9„ of the fourth-order equation ((55|) 
have the same sign, the four (complex or real) solutions 
of the real part is definitely negative. In this case, the 
general solution which is a linear combination of those 
four solutions, is stable. Since 94, 93, and 90 are positive 
defined, then the stability condition sought after is that 
92 and 9i must be simultaneously positive. The condition 
is 



1 



OR 



cq 



> 0. 



(61) 



"00=0 



Using thermodyanmical relation 



OPcq 



1 aPoc 

'T ds„ 



where 



,2 _ aPeq 



is the adiabatic velocity of sound, one 
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can rewrite the stability condition eq. (|6ip with the fonow- 
ing 



1 dP, 



oq 



T ds. 



1 

K 



(62) 



which is exactly the same condition found in previous work 
[21] in the Eckart frame. Note that, when k — ?> 0, one 
finds that the speed of sound can exceed unity violating 
causality. On the other hand, when k is a finite number 
restricted by 



1 dn 



cq 



1 1 aPcq 

< _ < 1 i '^'l 
K 



T ds. 



oq 



(63) 



then the velocity of sound satisfies < < 1. Thus, when 
the condition for k, eg. ((63)) is satisfied, the fiuid is stable 
against small perturbations and evolves without violating 
causality. 



4 Summary and concluding remarks 

We have proposed a novel formulation of the relativistic 
dissipative hydrodynamical model in an arbitrary frame 
by using extended matching conditions 

Uf.uJT'"' = A, Uf.SN'' = 5n. 

To apply the above extended matching conditions, we have 
also generalized the form of the entropy current for non- 
equilibrium state [cf. eq. (PD|) ] : 



The phenomenological parameter x introduced in the gen- 
eralization of the entropy current can be fixed by the ex- 
tended thermodynamic stability condition, 

/ib(5n + A~ xn = 0. 

(Note that in the usual formulation x = 0, because of 
A = and Sn — 0). Taking the thcrmodynamical stability 
condition into account, the entropy current is given by 



As seen in the above equation, the last term xnu^^ is 
the new correction term. The corresponding entropy pro- 
duction evidently does not depend on the Lorentz frame 
considered. It is given by 

In this paper, we chose the Landau frame and considered 
a dissipative fluid with zero chemical potential = 0, for 
simplicity. For this case, the % was assumed as 



x = K + ^n + e- 



The physical meaning of k, ^ and ^" are revealed in the 
discussion on the stability and causality of the fluid in 
Sec|31 The k is related to the bound for the speed of sound 
wave of the fluid and ^ and are related to the relax- 
ation of the off-equilibrium system. In the linearized field 
equations, the speed of sound Cg is actually restricted so 
that < Cg < 1 when k is chosen in the following range: 



1 dR 



cq 



T ds. 



cq 



1 

< - < 1 



1 dP, 



oq 



T 



cq 



On the other hand, the relaxation time of the perturbed 
n and ■K^'^ fields are respectively given by 

Tn = (2e + A^/i7o)C, = 2C"?7. 

The parameter k, also contributes to the relaxation time 
and may bring in a contribution to the initial condition 
characterized by TTq . 

We therefore conclude that, when the matching condi- 
tions 



A = xn = nn + ^n^ + f v 



and 5n = Q (also = 0) are imposed, the relativistic 
dissipative fluid can be applied to an analysis of the phe- 
nomena observed in the ultra-relativistic heavy-ion colli- 
sions not only in the Eckart frame, as discussed in |21) . 
but also in the Landau frame. The conditions that should 
be imposed seem to be independent of the Lorentz frame 
used, i.e.. 



AT 



and c„ < 



1 



1 9Po, 



oq 



(64) 



n 



where A and are thermal energy conductivity and relax- 
ation time of the thermal energy conduction, respectively 
|21| . In the Landau frame, should be regarded as 
because of the deflnition of the frame, = 0. 



We gratefully acknowledge discussions with T. Koide. The au- 
thor would like to warmly thank dr E. Infeld for reading this 
manuscript. 
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